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The XMM-Newton observation of GRB 040106 : evidence for an 

afterglow in a wind environment. 
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Abstract. We present the XMM-Newton observation of GRB 040106. From the X-ray spectral index and temporal decay, 
we argue that the afterglow is consistent with a fireball expanding in a wind environment. A constant density environment is 
excluded by the data. This is one of the very few cases in which this conclusion can be drawn. 
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Since the BeppoSAX revolution on the Gamma-Ray Burst 
(GRB) studies and the discovery of the GRB afterglows 
( Cost a et al. Il997l) . a canonical model has emerged to explain 
the afterglow properties : the fireball model dRees & Meszaros I 
Il992l iMeszaros & Rees Ill997l IPanaitescu et al. Ill998l) . This 
model is based on a blast wave which propagates into a sur- 
rounding medium, first consider ed to be uniform (a review 
of this case is done by IPiran I fl999). But this model (re- 
ferred as the InterStellar Medium [ISM] model) was not able 
to explain all the features in the afterglow spectra and light 
curves. First, some afterglow light curves displayed an achro- 
matic break (e.g. IPianetal.lEooi . This was interpreted by 
the non isotropy of the bla st wave , and t his refined mode l 
was called the Jet model ( lRhoadslll997t ISari et al. l fl999). 
Second, the optical afterglow light curves sho wed in some cas e 
a bump, associated with type Ic supernova (Reichart 1999). 
These and X-ray features (e.g. iReeves et al. Il2002t IPiro et al. I 
1 19991) show that long G RBs may be link ed with hypernovae 
and star forming region dMeszaros 11200 ll) . The density of the 
surrounding medium then decreases with the square of the dis- 
tance to the cent ral engine, due to the wind arising from the 
GRB progenitor (Cheva]iex^Xij2QQ0|) . This model is referred 
as the Wind model JDai & Lu I Il998t IMeszaros et ain ri998: 
IPanaitescu et al. ll998HChevalier & Li Il999l) . 

The use of the multi-wavelength observations of a GRB af- 
terglow allows one to determine the model parameters and to 
indicate the best model for each burst (and possibly to indicate 
why some bursts are best fitted with an ISM model while oth- 
ers are favoring a wind model). In this Letter, we will do so for 
GRB 040106, and we will show that the afterglow data of this 
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burst indicate a fireball interacting in a wind. We will discuss 
this fact and the implications. 



2. GRB040106 

GRB040106 was detected with the INTEG RAL Burst Alert 
Syste m on 2004 January 6 at 17:55:12 UTC llMereghetti et al. I 
2004). It was detected in the IBIS/ISGRI data as a ~ 60 sec- 
onds event. The peak fl ux was 6.5 X 10~ 8 erg ctrT 2 s _1 in the 
20-200 keV band JGotz et al. l2004 . A 45 kiloseconds long ob- 
servation with XMM-Newton began at 23:1 1:23 UTC on 2004 
January 6. The EPIC instruments were operating in full frame 
mode, with THIN filters (PN and MOS2) or MEDIUM filter 
(MOS1). One fading source was dete cted with in the error box 
of the INTEGRAL detection (Ehle et al. 2004), and was asso- 
ciated with the X-ray afterglow of GRB040106. Using a cross- 
correlation with USNO-A2.0 stars, the refined position of this 
source is RA: 1 1 52 1 2.43, Dec: -4647 15.9 ( J2000.0, 1 cr total 
positional error 0.7", iTedds & Watson I2004I) . 

The galactic position of this source is 1 = 292.5 and b 
= 14.88. In that direction, the column of density is N# = 
8.6 x 10 20 ctrr 2 dPickev & Lockmanlll990h. The gal actic op- 
tical extinction is E(B-V) = 0.1 dSchlegel et al. Il998l) . Optical 
observations detected an afterglow to this GRB, with an R mag- 
nitude of 22.4 + 0.1 on 2004 January 7 at 08:33 and 23.7 ± 0.3 
on 2004 January 8 at 08:25 (Masetti et al. 2004). This imply a 
R magnitude corrected for the absorption of our galaxy of 22. 1 
+ 0.2 and 23.4 + 0.4 (we assume a conservative 0.1 mag error 
in the reddening value). 

A radio observation detected a s ource not consistent with 
the position of the X-r ay afterglow dWieringa& Frail I [2004: 
ITedds & Watson 1 120041) . A later radio observation (on 2004 
January 21.48 UT) with the VLA gave two upper limits of 100 
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Fig. 1. PN light curve of GRB040106 in the 0.2-12.0 keV band. 
The light curve is shown with the best fit power-law model. The 
flux indicated at the right of this figure is given in the 2.0-10.0 
keV band. 

yuJy an d 170 aJv at 8.46 and 4.86 GHz respectively (2 <x upper 
limits, iFrail et al. l2004 . 

3. Data analysis 

We retrieved the Observation Data Files on the web page of 
the XMM-SOC, and processed them with the SAS 6.0 and the 
latest calibration files available. We have reduced the data us- 
ing the epchain and emchain task of the SAS. Due to high 
background activity, we discarded 8 kiloseconds of observa- 
tion, and kept only 37 kiloseconds of data. We have finally fil- 
tered the event files for good events (FLAG —— in the PN 
files, #XMMEA_EM for the MOS ones) and single or double 
patterns. We have extracted light curves and spectra using these 
filtered files, with a circular extraction region, with radius of 
45 arcseconds (MOS) or 25 arcsecond (PN, avoiding thus any 
CCD gap). Background light curves and spectra were extracted 
using a larger (4 times the surface of the source extraction re- 
gion) circular region on the same CCD. 

Because this afterglow is bright (more than 10 000 photon 
detected), we decided to use the chi square statistic and binned 
the spectra in order to obtain at least 20 photons from the source 
(net photons) in each bin. We present the light curve in Fig.^ 

Fitting the light curves with a canonical power law, we ob- 
tain a decay index of 6 = 1.4 ± 0.1 (we assume the notation 
F y oc v~ a t~ s , all errors are quoted at the 90 % confidence level). 
Such a decay is very common for a GRB afterglow. 

We then made a spectral study of this afterglow. We first 
used an absorbed power law model on the whole exposure (pre- 
sented in Fig.EJ, fitting both EPIC-MOS and EPIC-PN spectra 
simultaneously, and allowing a normalization to vary between 
each instrument to take into account cross-calibration uncer- 
tainties. The best fit values are indicated in Tabled We also 
looked for spectral variations. We divided the exposure in 3 
segments with different durations (6, 12 and 19 kiloseconds of 
live time respectively), and repeated the spectral study. All the 
results are displayed in Tabled 



Fig. 2. EPIC spectrum of GRB040106. The spectrum is shown 
with the best fit power law model. 



As one can see, the spectral parameter errors at the 90 % 
confidence level are compatible with no spectral variation. 

In this paper, we are primarily interested in the properties 
of the continuum. We thus carried out a simple analysis for 
line detection. We added to the power law model an individual 
gaussian line with a narrow width and free energy. We searched 
in the summed and the three slices of pn spectra. We did not 
find strong statistical evidence of line features. The most sig- 
nificant features are found at 0.64 keV (/ = (10 + 3) x 10~ 6 ph 
cm" 2 s- 1 ), 0.86 keV (/ = (5 ± 1.5) x 10" 6 ph cirr 2 s" 1 ) and 3.6 
keV (/ = (3.8 + 1.3) x 10~ 6 ph cirr 2 s _1 ) with an confidence 
level of 99.6%, 99.2% and 98.7% respectively, as derived from 
the F-test. 1 



4. Constraints on the fireball model parameters 

Using both the spectral and temporal properties of this af- 
terglow, we can investigate the parameters of the fireball 
model. We first tried to discriminate the type of fireball model 
(Wi nd, ISM or Je t m odels), using the closur e relationships of 
ISari et al. I dl998l) and lChevalier & Li I (1 19991) . We derived val- 
ues of 5 - 1.5a = 0.64 ± 0.13 and 5 - 2.0a = 0.42 + 0.14 (90 
% confidence level). The value of S - 1.5a should be used for 
the ISM and Wind models, while the value of 6 - 2.0a applies 
to the Jet case. We should expect values of -0.5 or for the 
the ISM model and -0.5 or 0.5 for the wind model, depending 
on the location of the cooling frequency, v e . We thus conclude 
that our values are compatible only with a wind model. In the 
jet model, we should expect values of or 1, which are not 
compatible with our values. 

The closure relationship for the wind model are 6 - 1 .5a = 
-0.5 if v c < Vjt and 6 - 1.5a = 0.5 if v c > v x . Thus, we got a 
compatibility only if the cooling frequency is above the X-ray 
band. 



1 a more precise assessment of the statistical significance as indi- 
cated by Protassov et al. (2002) would go beyond the scope of this 
paper 
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Table 1. Spectral parameters of GRB040106. Each segment is indicated with the corresponding time range (in second) from the 
burst in the observer frame 





Spectrum 


N H (10 20 cm" 2 ) 


Energy spectral index 


Flux (10 12 erg s 1 cm 2 ) 




d.o.f. 




Complete 


8.8 ±0.8 


0.49 ± 0.04 


1.01 ± 0.02 


1.07 


515 


1* 


segment (21185-27178) 


10 ±3 


0.5 ±0.1 


1.7 ±0.1 


1.06 


64 


2 nd 


segment (27178-41512) 


8±2 


0.43 ± 0.05 


1.25 ±0.02 


1.16 


271 


3 rd 


segment (41512-63019) 


8±2 


0.5 ± 0.05 


0.64 ± 0.02 


0.99 


227 



May the cooling frequency pass through the X-ray band 
during the observation ? The light curve should show a steep- 
ening at that moment, with a decay variation of 0.25. There is 
a deviation in the light curve at about 12 kiloseconds (in the 
XMM-Newton observation time), which degrades the fit (xl = 
2.37). We have tried to fit a broken power law to this light curve. 
The fit is not improved, indicating that the deviations are likely 
due to shot-time scale variations rather than to an overall steep- 
ening in the power law decay index. Moreover, the spectral in- 
dex should indicate a flattening with a variation of 0.5 when 
the cooling frequency passes through the X-ray band. We can 
rule out such a large variation of the spectral index. Finally, our 
spectral index and temporal decay values are compatible with 
a cooling frequency above the X-ray band in all the three seg- 
ments. We thus conclude that we indeed observe an afterglow 
in the wind model with the cooling frequency above the X-ray 
band. 

Using the theoretical model, we can now constrain some 
model parameters. The theoretical temporal and spectral slopes 
are 6 — {3p - l)/4 and a — (p - l)/2 respectively. We obtain 
p = 2.2 + 0.2 and p = 2.0 + 0.1 respectively. 

We can also use the time at which the cooling fre- 
quency passes through the X-ray band (using an upper limit) 
to constrain the surro unding density medium. According to 
IChevalier & Li I (1200 01 the cooling time observed in the X-ray 
band is : 

f ^ 1L5xl0 ^i±i) 3 (^-) 3 £- 2 ^Days (1) 

The redshift of this burst is unknown, we thus assumed the 
common median value of 1. Using the data from the prompt 
emission, we obtain E y ^2 ~ 1 ■ Because the cooling frequency 
is above the X-ray even in the first part of the observation, t c < 
0.23 days. We thus obtain : 

A*<2.0xlO- 3 (^J 3/4 £^ 4 (2) 

We now use the broadband spectrum between the radio and 
the X-ray band. The optical decay is 1.2 + 0.4, fully consistent 
with the X-ray decay. Also, the unabsorbed optical-to-X-ray 
spectrum is compatible with a single power law, independently 
supporting the conclusion that v c is above the X-ray band. In 
addition, this indicates that v m is below the optical band at the 
date of the first optical obser vation (0.61 days) . We use the ex- 
pression of v m given by Chev alier & Li I d200C)h : 

v m ~1.0xlO-(^) 2 (^) 1/2 ^Hz (3) 




1 o~ 6 1 , i , i , i , i , i I , 

108 ,Q10 10 12 , 14 10 18 10' 8 10™ 

Frequency (Hz) 

Fig. 3. Spectrum of GRB040106. The spectrum is shown at 
14.12 days after the burst. The power law indicated is the best 
fitted power law model from the X-ray data. 

Imposing v m < 4.28 x 10 14 Hz, we obtain : 

e^O.l^E^ (4) 

This condition does not give a very strong constraint on the 
value of e e : if eg is equal to one, then e e is less than ~ 0.4. 

Another constraint can be set by the flux density value in 
the X- ray. Using the equations given in Panait escu & Kurriarl 
(2000), we got at 3 keV (where absorption is negligible) and at 
the time of the XMM-Newton observation : 

'-"Met) (et) £ °™" ,y <5) 

The flux density measured by XMM-Newton is 6.56 x 
10~ 2 /Jy. We thus obtain : 

We have also verified if the position of v m can be con- 
strained by using radio data. We have extrapolated the X-ray 
and optical fluxes at the date of the radio observation, assum- 
ing that there is no jet, and thus no achromatic break in the 
light curves. As one can see in Fig|3] the radio upper limits are 
compatible with the optical/X-ray spectrum. The pres ence of a 
jet w ould produce a steepening in the light curves jRhoads I 
1997), and thus lower optical/X-ray fluxes. This would also 
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give a compatibility between the optical/X-ray fluxes and ra- 
dio upper limits, preventing us to derive other constraints on 

5. discussion and conclusion 

From the X-ray observations of the afterglow of GRB 040106, 
we have constrained the surrounding medium of the fire- 
ball to be a wind environment. We have derived a constraint 
on the density of the medium, which should be A, < 2 x 
10- 3 (e B /0.1)-°- 75 £^ 4 . By definition, if there is no radiative 
losses, then £52 = Ey^l^y, where e y is the conversion effi- 
ciency factor of the energy of the fireball in to gamma ra ys. We 
adopt here the value of e r ~ 0.2 as in Fr ail et al. 

Several GRB afterg lows have been observ ed to be possi- 
bly in the wind model ( Chevalie r et al. 1120 04). although only 
in a very few cases the wind profile is the only acceptable 
model (e.g. GRB01 1 121). In most of the cases, the A* value de- 



rived from the observations is ~ 0.3 - 0.7, compatible with the 
Wolf Rayet star wind observed in the Galaxv JChevalier & Li 1 
1 19991 IChevalieretal.112004. But the cases of GRB 020405 


(A. < 0.07|Chevalieret al. 


2004), GRB 021211 (A, = 0.0005 


Kumar & Panaitescu 2003; 


IChevalieretal.112004 and GRB 


011121 (A, = 0.003 IPrice et al. 1 


2002: Piroetal.||2004|) in- 



dicate a value of A, significantly lower than one. From present 
data, and using spectral modeling, we can not significantly con- 
strain A,. If we set e B = 0.1, e e = 0.3, then A, = 2.8 x 
10~ 3 . On the other hand, v alues of eg = 10~ 4 (observed by 
Panaitescu & Kumar 1l2002l in some afterglows) and e e = 0.3 

im ply that A» = 0.53 . 

IChevalier et all d2004 propose several interpretation for 
the low density sometime observed : a lower metalicity or a 
lower mass of the Wolf Rayet progenitors of GRBs, a GRB oc- 
curring along the progenitor rotation axis, and an unusual pop- 
ulation of Wolf Rayet stars responsible for some GRBs. This 
could be tested by observations of a large set of Wolf Rayet 
stars and the host galaxies of GRB afterglows surrounded by a 
wind. Other XMM-Newton afterglow observations could grow 
the sample of GRBs surrounded by a wind with a good known 
location, and thus allow one to test these hypothesis. 
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